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i-^ ■ Abstract 

We investigated the evolution of interacting disk galaxies using high-resolution iV-body/SPH simula- 
tions, taking into account the multiphase nature of the interstellar medium (ISM). In our high-resolution 
simulations, a large-scale starburst occurred naturally at the collision interface between two gas disks at 
the first encounter, resulting in the formation of star clusters. This is consistent with observations of 
interacting galaxies. The probability distribution function (PDF) of gas density showed clear change dur- 
ing the galaxy-galaxy encounter. The compression of gas at the collision interface between the gas disks 
first appears as an excess at ?ih ~ 10 cm~ 3 in the PDF, and then the excess moves to higher densities 
[~^. ' (tt-h ^ 100 cm~ 3 ) in a few times 10 7 years where starburst takes place. After the starburst, the PDF goes 

\Q , back to the quasi-steady state. These results give a simple picture of starburst phenomena in galaxy-galaxy 

t-H ■ encounters. 

Key words: galaxies:starburst — galaxies:ISM — ISM:structure — method:numerical 

in ■ 

o : 

OO , 1. Introduction low 10 4 K. (see §2 in Saitoh et al. 2008; hereafter pa- 
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per I). In most of early simulations, an isothermal equa- 
Galaxy-galaxy interactions (collisions and mergers) are tion of state (EOS) of 10 4 K was used for the ISM (e.g., 
believed to play important roles in shaping present-day Mihos & Hernquist 1996). Even when the energy equation 
galaxies, since they occur frequently in the standard hi- with a cooling term was used, the minimum temperature 



■ erarchical model of galaxy formation (e.g., White & Rccs was set to be ~ 10 4 K (e.g., Barnes & Hernquist 1996). 
1978; Lacey & Cole 1993; Kauffmann & Chariot 1998; Some authors (Springel 2000; Cox et al. 2004) adopted 
Okamoto & Nagashima 2001). Interactions drive not the EOS that became harder in higher density. They ar- 
only morphological transformations of galaxies but also gued that such EOS mimicked the effect of turbulence mo- 
bursts of star formation (e.g., Young et al. 1986; Barnes tions caused by SNe feedback from massive stars. Struck 
& Hernquist 1992; Smith ct al. 2007; Li ct al. 2008). (1997) employed a multiphase ISM model involving a low- 
A number of researchers have studied roles of galaxy- temperature phase under 10 K, although the modeling of 
galaxy interactions using numerical simulations (e.g., ISM was rather simple and mass resolution was low. 
Toomre & Toomrc 1972; Hernquist 1989; Mihos & Low-resolution simulations did not work well un- 
Hcrnquist 1994; Mihos & Hernquist 1996; Barnes & der the situation of strong galaxy-galaxy interactions. 
Hernquist 1996; Springel 2000; Cox et al. 2004; Springel Observations of interacting galaxies often show shock- 
& Hernquist 2005; Cox et al. 2006). Many simulations in- induced, widespread star-formation activities involving 
corporated the interstellar medium (ISM) and models for star cluster formation (e.g., the Antennae and the Mice 
star formation and supernova (SN) feedback from mas- galaxies) , whereas previous simulations showed that most 
sive stars. In such simulations, the numerical resolution is of star formation activities was concentrated to the central 
rather limited. The typical gas mass resolution is about regions of the galaxies. Barnes (2004) demonstrated that 
10 6 Mq. Because of this limited mass resolution, it was the collision simulations with the isothermal ISM and the 
impossible to treat the formation of cold gas phase be- star-formation model governed by the local energy dissi- 
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pation rate in shocks can reproduce large-scale starburst 
formed at the collision interface during the first collision. 
However, his shock-induced star formation model intro- 
duced an additional parameter and required the fine tun- 
ing of model parameters for star formation. 

The Schmidt law (Schmidt 1959) has been used as star 
formation model; 

dp* _ Q Pgas /j\ 
dt tdyn 

where p* and p gaB are the stellar density and gas den- 
sity, respectively, tdyn = (47rG , p gas )~ 1 / 2 is the local dy- 
namical time, and C* is the dimensionless star formation 
efficiency Since the isothermal EOS prevents ISM to form 
small-scale structures, in most simulations, stars had to be 
formed in fairly low-density gas. Recent models employ 
the threshold density for star formation, ~ 0.1 cm" 3 , 
above which star formation occurs. SNe releases energy 
into the surrounding ISM in the form of thermal or kinetic 
energy, or internal turbulence motions. When the value of 
C* was appropriately tuned, simulations based on these 
models successfully reproduced observed global properties 
of star formation in local disk galaxies, in particular the 
Schmidt-Kennicutt relation (Kennicutt 1998). However, 
simulations with the isothermal EOS or the temperature 
cutoff in cooling functions at 10 4 K failed to reproduce 
three-dimensional structure of star formation activities, 
simply because the star formation in these simulation took 
place effectively everywhere in the gas disk. In paper I, wc 
demonstrated that high-resolution simulation of the ISM 
with low temperature cooling below 10 4 K, combined with 
higher threshold density for star formation, gave rise to 
thin and inhomogeneous star forming regions as observed 
in isolated disk galaxies. Our high-resolution simulations 
gave a good description of small-scale structures of star 
forming regions, and at the same time well reproduced 
global characteristics such as the Schmidt-Kennicutt rela- 
tion. 

In this paper, we report the results of high-resolution 
simulations of a collision between two disk galaxies with 
realistic ISM with cooling term and high-density star- 
formation criterion for the Schmidt-type star formation 
model. This ISM model is based on the one used in pa- 
per I and it can reproduce properties of star-formation in 
local disk galaxies, without the fine tuning of the value 
of C*. Our main conclusion is that high-resolution sim- 
ulations can resolve shock-compressed regions, and, as a 
result, these simulations can naturally reproduce a large- 
scale starburst at the collision interface of gas disks during 
the first encounter. Initially, the density of the gas in the 
shock-compressed region reaches nn ~ 10 cm -3 . This com- 
pressed gas collapses through radiative cooling and gravi- 
tational instability in the timescale of ~ 10 7 yr and reaches 
a high density at which star formation is triggered. Thus, 
our modeling of the ISM, which requires sufficiently high- 
resolution, naturally reproduces the shock-induced star- 
burst. A comparison run, which was performed with the 
temperature cutoff in cooling function at 10 4 K and rel- 
atively low-density threshold for star-formation, displays 



only a small enhancement of the star formation rate (SFR) 
during the first encounter. By using high-resolution simu- 
lations with a realistic star formation model, we can suc- 
cessfully reproduce not only the quiescent star formation 
in disk galaxies (paper I) but also the large-scale starburst 
at a collision interface and nuclear starbursts observed in 
interacting galaxies without any extra assumptions. 

2. Method and setup 

We used a parallel tree SPH code ASURA (Saitoh, 
in preparation) that utilizes the special-purpose hard- 
ware GRAPE. We first constructed near-equilibrium, self- 
consistent galaxy models with a dark-matter halo and an 
exponential disk, and then we converted a part of parti- 
cles into gas particles. The TV-body realization was gen- 
erated by GalactlCS (Kuijken & Dubinski 1995). Our 
galaxy model is a low-mass, gas rich Sc galaxy with the 
characteristic rotation velocity of 120 km s _1 . It has a 
gradually increasing rotation curve, which is consistent 
with observations of nearby low-mass galaxies (Sofue et al. 
1999). The dark-matter halo mass and the disk mass are 
1.05 x 10 11 M Q and 6.3 x 10 9 M 0! respectively. The scale 
length of the exponential disk is set to be 4 kpc (Roberts 
& Hayncs 1994). We converted 20 % of the disk particles, 
1.2 x 10 9 Mq, into gas particles with the exponential pro- 
file of an 8 kpc scale length. We also converted 1 % of 
the dark matter particles into gas particles with the same 
profile as that of the dark matter halo. The initial gas 
temperature is set to be 10 4 K. The galaxy-galaxy col- 
lision is characterized by a prograde parabolic orbit with 
the pericentric distance i?p Cr i = 7.5 kpc and the initial sep- 
aration of 75 kpc. The relative velocity at the moment of 
the pericenter passage is ~ 370 km s _1 for the point-mass 
approximation. Two galaxies have the same spin direc- 
tion. We prepare three initial models with different mass 
resolutions. Further details of the set up are described in 
Saitoh et al. (in preparation). The gravitational soften- 
ing was set to be 20 pc for all particles in all runs. The 
particle mass and number of particles are summarized in 
table 1. 

We use a cooling function for a gas with the metallicity 
of half the solar value for the temperature range from 10 K 
to 10 8 K (Spaans & Norman 1997; Wada & Norman 2001). 
A uniform heating from the far-ultraviolet radiation is in- 
cluded. The intensity of the far-ultraviolet radiation is 
set to be the half of the solar neighborhood value (Wolfirc 
et al. 1995). 

Models for the star formation and the SNe feedback are 
similar to those in paper I. We first adopt the model which 
allows a wide range of temperature for ISM (10 — 10 8 K) 
and uses the values 100 cm -3 and 100 K as the threshold 
density and temperature for the star-formation, respec- 
tively. We add the label 'C in run names for these models. 
For comparison, wc perform a run which only allows the 
ISM temperature above 10 4 K and employs 0.1 cm" 3 and 
15000 K as the threshold density and temperature for the 
star-formation, respectively The behavior of the ISM in 
this run should be similar to those in previous simulations 
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with isothermal ISM of ~ 10 4 K. The label 'A' indicates 
this model. We adopt C* = 0.033 for all runs. 

3. Starburst during the first encounter 

Figure 1 shows global SFR during the first encounter. 
Before the encounter, each galaxy has a nearly constant 
SFR of - 1 M Q yr" 1 . For Runs MIC, M2C, and M3C, 
SFR increases by almost an order of magnitude just after 
the pericenter passage (~ 430 Myr; see the thin curve). We 
call this peak as a starburst. This duration of the starburst 
is short and the phase of star formation returns quickly 
to the quiescent star-formation phase after the first pas- 
sage. The slightly larger SFR than those seen in the pre- 
collision might be explained by strong non-axisymmertic 
perturbations from the companion galaxy (Noguchi 1988). 
The behavior of the starburst at the first encounter is sim- 
ilar to that obtained by Barnes (2004) who implemented 
a shock-induced star formation. Also, the difference be- 
tween Runs MIC, M2C, and M3C is very small. We con- 
clude that our result is independent of the mass resolution 
as long as the resolution is sufficiently high. In contrast to 
these 'C runs, Run Ml A shows very small enhancement 
of SFR during the first encounter. 

The difference in SFR is owing to the hydrodynamic 
pressure of the gas that forms a bulk of the giant filament. 
The pressure of gas in Run MIA is roughly a thousand 
times larger than those in Runs MIC, M2C, and M3C 
at the same density, because the gas in Run MIA can 
not cool below 10 4 K while those in Runs MIC, M2C, 
and M3C can cool down to 10 K. In Run MIA, the high 
pressure ISM in the contact interface resists collapse even 
under significant compression. As a result, the density 
of the ISM is not increased significantly at the interface. 
Since the timescale of the pericenter passage of the galaxy- 
galaxy interaction is not very long (say, 10 7 years), in Run 
MIA, a smaller amount of the stars are formed during the 
first encounter, as long as we adopt the Schmidt law for 
the modeling of the star formation. 

Figure 2 shows the gas density of Run M2C at the or- 
bital plane for six different epochs. Two galaxies first fol- 
low the parabolic orbit and then collide at the pericenter 
(top panels). During the first encounter, the maximum 
Mach number reaches ~ 100 for the cold (T < 100 K) gas. 
The ISM at the collision interface is strongly compressed 
by shocks and forms a large and dense filament (see the 
top-right panel of figure 2; t — 420 Myr). We call it a giant 
filament. The length of this giant filament at its formation 
epoch is ~ 10 kpc. We measured the gas mass in the rect- 
angular box of the size 10 kpc x 2 kpc x 2 kpc along the 
giant filament (the black rectangular box in the top-right 
panel of figure 2). The gas mass of the giant filament is 
~9x 10 s Mq. This mass is about one quarter of the total 
gas mass in the system. The stellar mass formed during 
the starburst phase (t = 420 - 450 Myr) becomes ~ 20 % 
of the initial filament mass. The distribution of star form- 
ing regions is consistent with observations in interacting 
galaxies (e.g., Whitmore & Schweizer 1995). Overall evo- 
lution of the colliding galaxies in the first encounter is 
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Fig. 1. Global SFR during the first encounter as a func- 
tion of time. Thick solid, dotted, dot-dashed, and 
dashed curves indicate the global SFR for Runs MIC, 
M2C, M3C, and MIA, respectively. The thin curve 
indicates the separation between two galactic centers. 



quite different from those in previous simulations with the 
isothermal ISM and the Schmidt-type star formation mod- 
els. Our high-resolution simulation naturally captures the 
shock-induced, large-scale starburst at the collision inter- 
face. A number of "star clusters" form in the giant fila- 
ment at the starburst phase. We will study the formation 
of the star clusters in further details in forthcoming pa- 
pers. 

In addition to the shock-induced starburst in the gi- 
ant filament, we also find the starburst during a core- 
merging phase. Figure 3 follows the evolution of SFR 
during this phase. As is shown in this figure, simulations 
with our model shows starbursts with multiple peaks in 
the final core-merge phase. The star formation mainly 
takes place within compact clumps in this phase. This 
is consistent with the observations of ULIRGs (Sanders & 
Mirabel 1996). Our simulations successfully reproduce the 
nuclear starbursts as well as the large-scale starbursts. We 
further find the star formation enhancement in the tidal 
tails. 

4. Probability distribution function during the 
starburst 

Figure 4 shows the volume weighted probability dis- 
tribution function (PDF) of Run M2C for several differ- 
ent epochs before, during, and after the first encounter. 
Although there are some noticeable changes in the PDF 
with time, they are difficult to see because of the large 
vertical scale. In order to see the changes in the den- 
sity distribution of the gas more clearly, we plot the mass 
functions of the gas, dM / dlog(nn) , in figure 5. Just be- 
fore the encounter (t = 410 Myr), an excess of the gas 
(bump) appears at nn ~ 10 cm" 3 . By t — 430 Myr, this 
bump moved to tih ^ 100 cm" 3 , and by t = 450 Myr it al- 
most vanished. Note that this period of t = 430 — 450 Myr 
coincides with the period of the starburst. Thus, we 
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Table 1. Number and mass of particles for initial conditions of merger simulations. The cutoff temper- 
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for star 




formation arc 
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Run 
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N Gas ° 
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rp c 
J cut 
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MIA 


6 930 000 


341 896 
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M Q 
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0.1 cm" 


-3 


15000 K 


0.033 


MIC 


6 930 000 
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3 x 10 4 


M 


10 K 
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-3 


100 K 


0.033 


M2C 


13 860 000 
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296 322 
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M3C 


27 720 000 
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a The number of DM particles. b Thc number of disk particles. c The initial number of SPH particles. d Mass of 
individual particles (M ). c Cutoff temperature of cooling function (K). Threshold density of star formation (cm -3 ). 
s Threshold temperature of star formation (K). h Star formation efficiency. 




Fig. 2. Gas density distribution for Run M2C. Each panel shows a 16 kpc X 16 kpc re- 
gion. The elapsed time in Myr is shown the right-bottom corner of each panel. 



have a simple and clear picture of the shock-induced star- 
burst; Initially, large amounts of the gas is compressed to 
nn ~ 10 cm -3 . Then, this compressed gas cools and be- 
come denser through gravitational instability, and when 
the gas becomes dense enough to form stars, the star- 
burst takes place. The starburst is led mainly by star for- 
mation in massive clusters. After the high-density gas is 
consumed, the density PDF returns to the quasi-steady 
state. In figure 6, we show the time variation of the 
mass in selected density ranges. The relative changes 
are larger for higher densities. Moreover, the peak is 
found at systematically later epochs for higher densi- 
ties. The time delay from = 10 cm~ 3 to 100 cm~ 3 
is 10 Myr. This timescale corresponds to the local dy- 
namical time, £dyn(?T-H = 10 cm~ 3 ) ~ 10 Myr, and it is 



faster than the mean evolution timescale in isolated disk 
galaxies (~ 5 idyn! see figure 10 in paper I). The lifetime of 
the mass excess in each density range in figure 6 is around 
20 Myr. As is expected, this timescale corresponds to that 
of the starburst. 

5. Summary 

We have studied the evolution of the multiphase ISM in 
the colliding galaxies. Our simulations, for the first time, 
resolved the shock-induced gas compression and starburst 
at the collision interface during the first encounter of col- 
liding galaxies, with the density dependent star formation 
model (the Schmidt law). The behavior of the ISM is 
qualitatively different from the results of previous studies 
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Fig. 3. Global SFR during the core-merger phase as a 
function of time. Thick solid and dotted curves indi- 
cate the global SFR for Runs MIC and MIA, respectively. 
For reference, the peak SFRs at the first encounters is 
12 — 16 Mq yr — 1 (the horizontal gray zone), which is small 
to the SFR peaks in the core-merger phase shown here. 
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Fig. 4. Volume weighted probability distribution func- 
tions for the gas of Run M2C for five different epochs. 
The gas within R < 10 kpc and \Z\ < 1 kpc, where 
R is the distance from the mid-point of two galax- 
ies at the initial time in the orbital plane and Z is 
the vertical height from the orbital plane, is considered. 



with isothermal ISM. 

Shock-induced starbursts, which have been observed in 
many interacting galaxies, naturally take place in our 
high-resolution simulations without any extra assump- 
tions for the star-formation, such as the dependence on 
the energy dissipation rate (Barnes 2004; Okamoto et al. 
2005). In previous simulations with limited resolutions, 
significant shock-induced starbursts took place only when 
the presence of the shock itself is used as the criterion for 
star formation. 

In paper I, we showed that high-resolution simulations 
with realistic ISM and high-density star-formation crite- 
rion can reproduce the Schmidt-Kcnnicutt relation with- 
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Fig. 6. Time evolution of gas mass for five different den- 
sity ranges. Long (green) and short (light blue) hori- 
zontal bars correspond to local dynamical times of den- 
sities at nn = 10 cm -3 and 100 cm~ 3 , respectively. 



out the need of the fine tuning of the value of C* . Previous 
simulations with isothermal ISM and low-density star- 
formation criterion required fine tuning. In this paper, we 
showed that high-resolution simulations can directly re- 
solve shock-compressed regions and these simulations can 
automatically reproduce shock-induced starbursts. High- 
resolution (particle mass less than 10 4 Mq) with a realis- 
tic star formation model allows us to study both quiescent 
star formation and starburst modes in iV-body/SPH sim- 
ulations of galaxy formation. 
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improve our manuscript. The author (TRS) thanks 
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